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Abstract - We investigate the combined effect of Hund's and spin-orbit (SO) coupling on su- 
perconductivity in multi-orbital systems. Hund's interaction leads to orbital-singlet spin-triplet 
superconductivity, where the Cooper pair wave function is antisymmetric under the exchange of 
two orbitals. We identify three (i-vectors describing even-parity orbital-singlet spin-triplet pairings 
among t2 g -orbitals, and find that the three d-vectors are mutually orthogonal to each other. SO 
coupling further assists pair formation, pins the orientation of the d- vector triad, and induces spin- 
singlet pairings with a relative phase difference of tt/2. In the band basis the pseudospin d- vectors 
are aligned along the z-axis and correspond to momentum-dependent inter- and intra-band pair- 
ings. We discuss quasiparticle dispersion, magnetic response, collective modes, and experimental 
consequences in light of the superconductor Sr2Ru04- 
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Introduction. — Since its inception, standard 
Bardeen- Cooper- Schrieffer (BCS) theory has been con- 
sidered a classic example for a collective phase emerging 
from quantum many body effects. However, the discovery 
of unconventional superconducting phases near antiferro- 
magnetic order in heavy fermion compounds [TJ[2] , organic 
materials [3], and, most recently, Fe-pnictides [3] have ex- 
posed the limits of a single-band BCS formulation. The 
origin and nature of superconductivity in complex mate- 
rials where multiple bands cross the Fermi level therefore 
remains a field of active research, harbouring intriguing 
challenges and mysteries. 

In particular, when the electronic structure near the 
Fermi energy is composed of different orbitals and spins 
mixed via spin-orbit (SO) coupling, a pairing symmetry 
analysis could be non-trivial. For example, a local micro- 
scopic interaction such as Hund's coupling may naturally 
favour inter-orbital spin- triplet pairing between electrons. 
However, when orbital and spin fluctuations are significant 
due to inter-orbital hopping and SO interaction, pairing 
in definite orbital and spin channels (e.g., spin-singlet or 
-triplet pairing between electron in orbitals a and b) is 
not well defined. Equivalently, from a Bloch band per- 
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spective, where the kinetic Hamiltonian including SO ef- 
fects is diagonal, the decoupling of the microscopic interac- 
tion effectively leads to intra- and inter-band pairing with 
pseudospin-singlet and/or -triplet character. 

Below we present a systematic study of how SO and 
Hund's couplings jointly give rise to superconductivity in 
t2 g (i.e., d yz , d xz , and d^) orbital systems. Our find- 
ings may apply to a number of multi-orbital d-subshell 
superconductors. To be specific we base our quantitative 
considerations on the proposed chiral spin-triplet super- 
conductor Sr2Ru04. Here, despite intense investigation 
for more than a decade, a clear picture for the pairing 
symmetry, the pairing mechanism and the relevant bands 
involved that is consistent with all experimental observa- 
tions has not yet emerged [5j[6]. 

The paper is organized as follows. In the second sec- 
tion we discuss Cooper pairing in multi-orbital systems. 
We find that superconductivity from local Hund's ex- 
change can naturally be characterized by three mutually 
orthogonal d-vectors each describing inter-orbital even- 
parity spin-triplet pairing. We then show how SO coupling 
pins the orientation of the d- vector triad and induces and 
enhances pairing via coupling to spin-singlet pairing order 
parameters with a fixed relative phase difference of it/ 2. 
In the third section, we map these local pairing order pa- 
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rameters, defined in an orbital and spin basis, to inter- 
and intra-band pairing in the Bloch band basis. Pairing 
in the Bloch bands has a strong momentum dependence 
and the magnitude and direction of the d- vectors depend 
on the orbital composition at each k-point. In the fourth 
section, we present the complete self-consistent mean-field 
(MF) results involving 9 complex order parameters using 
band structure parameters that reproduce the Fermi sur- 
face (FS) reported on Sr2Ru04. In addition, the resulting 
anisotropic quasiparticle (QP) dispersion, the magnetic re- 
sponse and the critical pairing strengths in the presence of 
SO coupling are considered. We summarize our findings 
and discuss the relevance for SO-coupled d-orbital super- 
conductors such as Sr2Ru04 in the last section. 

Pairing in SO coupled t2 g systems via Hund's in- 
teraction. — For multi-orbital 3d-subshell systems such 
as the Fe-pnictides, it was recognized that Hund's coupling 
(interaction strength denoted by J) is as important as on- 
site Coulomb repulsion (U) [HIE], while SO coupling (2A) 
is relatively weak [9]. In contrast, recent x-ray measure- 
ments on 5d transition metal compounds such as Ir-based 
oxide materials found that the SO interaction of 0.6 eV 
is roughly comparable to the on-site Coulomb energy [10] , 
suggesting that SO interaction is larger than Hund's ex- 
change (since J < U). Given that the effective pairing 
interaction in the spin-triplet channel arising from Hund's 
coupling and inter-orbital Hubbard repulsion (V = U—2J) 
scales as V — J = U — 3 J (see below), we therefore ex- 
pect that for 4d-subshell materials such as Sr2Ru04 both 
SO and spin-triplet pairing interactions are intermediate 
in strength and of similar magnitude [TT1IT7] . Since nei- 
ther interaction is negligible nor dominant, we treat both 
on an equal footing in the present study. 

While on-site Hund's and further neighbor exchange in- 
teractions have been recognized to be important for spin- 
triplet pairing [7, 18-21 j, the combined effect of SO and 
Hund's couplings on inter-orbital spin-triplet pairing has 
not been investigated in t2 g -orbital systems. To under- 
stand superconductivity in SO coupled t2 g -orbital sys- 
tems, we consider a generic Hamiltonian H = H^ n -\-Hso-\- 
H[ nt consisting of kinetic, SO, and local Kanamori inter- 
action terms. In this section we leave the kinetic Hamil- 
tonian i^kin unspecified and focus on the pairing proper- 
ties arising from the interplay of the atomic SO coupling 
Hso = 2A ~Li • Si and the local interaction, which, pro- 
jected on the t2 g orbitals, are given by 

H S o = iX^^eauc^c^a 1 ^,, (1) 

i abl 
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Here and in the following, summation over repeated spin 
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Fig. 1: (Color online) The orbital-singlet spin-triplet d- vectors 
form a triad whose orientation is pinned along x, y, and z (or 
— x, — y, and — z) in the presence of SO coupling. See main 
text for details. 

indices a, a' =t, I is implied while the indices a, b G 
{yz, xz, xy} belong to an ordered set of t2g-orbitals. Fur- 
thermore, a 1 stands for Pauli matrices, c^J creates an elec- 
tron on site i in orbital a with spin a, and e a u denotes 
the totally antisymmetric rank-3 tensor. For transparency 
we have also introduced separate interaction strengths 
for Hund's coupling (J) and pair hopping (J 7 ), although 
J = J' at the atomic level. 

Let us apply a MF approach to study the particle- 
particle instabilities of the microscopic interaction H[ nt 
using the following zero momentum pairing channels 

K a/b = ^E^]-'( C ^ C -^+ C k.C a _ k(T ,), (3) 
k 

= ^E[^^]-'( c ^ c -w- c k.c a _w),(4) 

k 

where N is the number of k points. Here, = 
(A^ 6 ) (= A^ a ) stands for intra- (a = b) and inter- 
orbital (a 7^ b) spin-singlet pairing, which is even un- 
der the exchange of orbital quantum numbers {i.e. they 
form "orbital triplets"). The vector order parameter 
d a /6 = ((d x a/b ), (d y a/b ): (d z a/b )) (= -d b/a ) on the other hand 
parametrizes inter-orbital (a ^ b) spin-triplet pairing con- 
sistent with the usual d- vector notation where z(d • a)a y 
describes the spin-triplet pairing gap |2,22j. Note that 
d a /5 is odd under orbital exchange, which is characteristic 
of an "orbital singlet" (while d a / a = 0). Note also that the 
above order parameters are all even under a parity trans- 
formation as they are locally defined; this feature differs in 
particular from conventional odd-parity spin-triplet pair- 
ing where orbital degrees of freedom are absent. 

Using the above pairing channels the interaction Hamil- 
tonian takes the form 

ffint KaK/a + (V ~ J)N £ /b 

a a,b,l 

+ J'N E K)a^t /b + (V + J)N E K% (5) 

a^b a^-b 

where it is clear that only Hund's coupling can give rise 
to an instability in a spin-triplet channel [7J[T9]. We thus 
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concentrate on the effective pairing interaction 
H! nt = (U-3J)Nj2d l l /b d l a/b 

a,b,l 



-.an 



(6) 



in the attractive regime U/3 < J (< U). In gen- 
eral, orbital-singlet spin-triplet pairing can also induce 
spin-singlet pairing so that the remaining terms in Eq. 
(0 would hamper spin-singlet pairing. However, we as- 
sume that their effect is negligible to keep the follow- 
ing self-consistent calculations feasible, and since the in- 
duced spin-singlet pairing amplitudes are for the most 
part smaller than the spin-triplet pairing amplitudes (see 
below). For notational clarity we label in the follow- 
ing inter-orbital pairing only by the three combinations 
a/b = xz/xy,yz/xy,yz/xz. 

To understand the effect of SO interaction, let us re- 
mark on pairing in the absence of SO coupling first. 
In the case of the layered compound considered below 
(and for a rather large parameter range) the three spin- 
triplet d- vectors d xz / xy , d yz / xy , and d yz j xz form a triad 
of mutually orthogonal vectors with an arbitrary orien- 
tation and chirality in spin space, and no relative com- 
plex phase difference (hence preserving time reversal sym- 
metry (TRS)). This can be understood by analyzing the 
Ginzburg-Landau (GL) free energy, which without SO 
coupling is given by 

T ~ X)[^l^l 2 + ^ 1) (d,-d:) 2 + B( 2 )|d I/ .d I/ | 2 ] 

+ ^[cW(d,.d I ,)(d K .d K r+c( 2 ; )|d I ,i 2 |d K | 2 (?) 
+ c% |d„ . d K | 2 + c< 4 » • d^i 2 + c^k . d:) 2 ] 

up to fourth order, by analogy to He-3 [23]. Here v,k 
stand for orbital pairs a/b, while the (real) quartic mix- 
ing parameters obey C^l — C^v and the asymmetry 
between in-plane and out-of-plane orbitals due to e.g. 
inter-orbital hopping is reflected in distinct coefficients 
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etc.). This form is dictated 



by gauge symmetry, SU(2) spin rotational symmetry, time 
reversal symmetry and the underlying lattice symmetries, 



and shows that the and C^l terms are sensitive to 
the relative orientation of the d- vectors, whereas the cjy^ 

(5) 

and Cvk contributions additionally depend on their rela- 
tive complex phases. 

However, once SO coupling is included, d xz / xy , 
d yz / xyi and d yz / xz are pinned along x, y, and z 
directions, respectively, as shown in fig. 1. In- 
version/time reversal symmetry on the other hand 
is still preserved and reflected in the degeneracy of 
the orientations/chiralities {d xz / xy , d yz j xy , d yz j xz \ and 
{-d xz/xy ,-d yz/xy ,-d yz/xz }. The pinning of the d- 
vectors occurs due to additional terms in the free en- 



ergy such as ~ a^\d 
b^\\d x , I s 

LI yz/xz\ 



yz/ J 2 + a (2) [K z/ J 2 

K /xz \ 2 ] +b^H z/x / + \di /xy 



d z l 2 l 

xz/xy\ J 



where 



] \d x (d y )* + d y (d x Y 

L yz/xy\ xz/xy' yz/xy^ xz/xy' 

the expansion parameters depend on the SO coupling 
strength, naively suggesting that < b^ 2 \ 

etc. \j SO interaction furthermore leads to a linear cou- 
pling between a particular component of (inter-orbital) 
spin-triplet pairing and (intra-orbital) spin-singlet pair- 
ing. For example, writing SO coupling between yz and xz 
orbitals in the form of — iA^]^/ (cj^. cj 



^xz _ xzryz x f , 

following linear coupling is allowed in the GL free energy: 
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Note that d yz j xz prefers the z-direction by coupling to 
spin-singlet pairing with a relative phase difference of 
±7r/2 depending on the sign of A. This is consistent 
with our findings below that the spin-triplet order param- 
eters are purely real while the spin-singlet amplitudes are 
purely imaginary. A similar analysis can be carried out 
for cP / „ and d y , . The overall order parameter for 

xz/xy yz/xy 

yz and xz orbitals then is d z xz/yz + i(A° xz/xz + A* yz/yz ). 
Since the relative phase between the orbital-triplet spin- 
singlet and the orbital-singlet spin-triplet order param- 
eters is fixed, there should be a collective mode rep- 
resenting a resonance of supercurrent flow between the 
coupled order parameters with an energy scale of order 
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Note that the above result is fundamentally different 
from similar two orbital models, which lead to a single 
orbital-singlet spin-triplet d- vector [19.21.24 j. The present 
model is also distinguished from other models where the 
momentum dependence in the band pairing usually orig- 
inates from nonlocal momentum dependent interactions 
[18], whereas here it arises from spin and orbital mixing 
in the Bloch bands as described next. 



Momentum- dependent pairing in the Bloch 
bands. — Despite having uniform pairing amplitudes 
d yz / xz , d yz/xy , d xz/xy , & s yz/yz , ... the corresponding inter- 
and intra-band pairings in the Bloch band basis (now car- 
rying band and pseudospin quantum numbers - 77, p = 
a,/?, 7 and s = ±) acquire a strong momentum depen- 
dence due to the mixing of orbitals through hopping and 
SO coupling. To understand how the above local pairing 
in the orbital and spin basis corresponds to pairing in the 
Bloch band basis, let us introduce the kinetic Hamiltonian. 
The most generic kinetic Hamiltonian for t2 g orbitals in a 



1 Analyzing the energetics of a corresponding two orbital model 
one can indeed show that SO interaction tends to stabilize e.g. the 
d z , -component over d x , or d y , . 

yz/xz r yz/xz yz/xz 
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Fig. 2: (Color online) Momentum- resolved pairing amplitudes 
in the Bloch band basis for 3 J — U — 0.9 and A = 0.15. Panel 
(a) and (b) represent inter- and intra-band pairing, respec- 
tively. The grey lines indicate the /?, 7, and a FS sheets (from 
inside to outside). Note that pairing from Hund's coupling 
preferentially involves electronic states near the FS sheets and 
that the intra-band pairing amplitudes are about one order of 
magnitude larger than inter-band pairing amplitudes. 



single layer perovskite structure has the form 
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where 

are ^l xz = — 2£iC0S& 2/ / x — 2t2Cosk x / y — /ii, £ k 
—2^3(008^ + cosky) — At4COsk x cosk y — H2, and ^ = 
— ^sin/c^sm/c^. For the MF calculation below we have 
chosen the parameters t\ = 0.5, £2 = 0.05, ts = 0.5, 
£4 = 0.2, £5 = 0.05, \i\ = 0.55, and \±2 = 0.65 (all en- 
ergies here and in the following are expressed in units of 
2t\ = 1.0). The underlying FS obtained from diagonaliz- 
ing i^kin with SO coupling strength A = 0.15 is shown in 
fig. [2] along with momentum-dependent band pairing am- 
plitudes. The FS agrees well with first principles calcula- 
tions [14J and the experimentally measured FS of S^RuC^ 
[T7II251 26 , consisting of three bands labelled a, /?, and 7. 

In the presence of SO coupling the bands are mix- 
tures of all three orbitals and different spins, e.g. £ k+ = 

ik c kf + ^k c kt + ^k c kl (V = a > At)- Hence consider- 
ing inter- and intra-band pairing amplitudes in the band 
basis, it is clear that the x- and ^/-components of the 
inter-band pseudospin-triplets such as (£k±£-k±) van i sn ? 
since (d£|<i^ k ^), (^kf^-k|)> and {d^d x _\^j amplitudes are 
zero (similarly for t^i)- Thus only finite z-components 
of the three inter-band pseudospin-triplet d-vectors and 
inter-band pseudospin-singlet order parameters (such as 
(£k+£-k- ^^k+^-k-)) can a PP ear - Figure [2] reveals that 
intra-band pairing is strongest and sharply peaked around 
the FS due to the mixing of all orbitals via SO interac- 
tion and inter-orbital hopping, and the ideal conditions for 
zero-momentum pairing. Inter-band pairing in contrast is 
about an order of magnitude weaker and, in particular for 
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Fig. 3: (Color online) MF solutions for different SO coupling 
strengths for the Sr2Ru04 based band structure. Orbital- 



singlet spin-triplet pairing d x 



'./xyi G-yz/xy: 



and d 



yz/x 



(purely 



real) induces finite intra-orbital spin-singlet pairing Ay Z / yz , 
^xz/xzi an d ^xy/xy (purely imaginary). We also checked for 
induced inter-orbital spin-singlet pairing amplitudes, which, 
however, vanish. 



(£k+£-k-)> more spread out in momentum space, marking 
Bloch band states that are energetically still close enough 
to the FS to participate significantly in pairing. 

This analysis demonstrates that inter-orbital pairing 
arising from Hund's interaction leads to k-dependent 
inter- and intra-band pairing in pseudospin-singlet and 
and pseudospin-triplet (z component only) channels. Fur- 
thermore, the pairing instability occurs simultaneously 
within and between all bands rather than in a single active 
band with superconductivity leaking into passive bands 
through, e.g., pair hopping. The role of intra-band spin- 
triplet pairing between a and /? bands in multi-orbital su- 
perconductors like Sr2Ru04 has also been the focus of re- 
cent studies, where the inter-band order parameter, how- 
ever, breaks TRS [27] and an intrinsic anomalous Hall 
effect can contribute significantly to a large TRS breaking 
signal in Kerr rotation experiments [28,29 . 

Pairing transition, QP dispersion, and magnetic 
response. — For concreteness we study the effect of SO 
coupling on spin-triplet pairing originating from Hund's 
interaction, including the QP dispersion and the magnetic 
response. As discussed in the previous sections the quali- 
tative results are generic for SO coupled t2g-bands (or p- 
orbital systems) and can be applied to specific materials 
such as the single layer ruthenate [5]l6] and the Fe-pnictides 
[7(130] using the appropriate band structure. 

Using the kinetic Hamiltonian of eq. (fTU|) with a pa- 
rameter choice mimicking the single layer ruthenate band 
structure, the MF solutions for various A are displayed in 
fig. [3j As one can see, in the absence of SO interaction 
an orbital-singlet spin-triplet pairing instability develops 
at a large coupling strength 3J — U > 1.0 for d xz / xy and 
dy Z /xy Although numerically difficult to resolve, we ex- 
pect that dy Z / xz and the intra-orbital spin-singlet order 
parameters simultaneously become finite through quartic 
or higher order couplings in the Landau free energy ex- 
pansion. While the magnitudes of the order parameters 
depend on the details of the band structure, a robust fea- 
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Fig. 4: (Color online) QP bands for 3J-U = 0.9 and A = 0.15. 
Panel (a) is a magnification of panel (c) about the Fermi level, 
revealing the gaps opening up on the FS sheets. Panel (b) 
shows the DOS and the QP gap near the Fermi level. 
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Fig. 5: (Color online) Magnetization parallel to the applied 
magnetic field B for A = 0.075 (top) and 0.15 (bottom) and two 
field orientations at 3 J — U = 0.9. The solid lines represent to- 
tal magnetization, dashed lines stand for orbital contribution, 
and dash-dotted lines for spin magnetization. For sake of com- 
parison the magnetic response both in the presence (orange) 
and in the absence (grey) of superconductivity is displayed. (B 
is expressed in units of 2ti = 1.) 



ture is that finite SO coupling drastically reduces the crit- 
ical pairing strength. This reduction is mostly facilitated 
by the additional hybridization provided by i^so, which 
helps to overcome the momentum mismatch between or- 
bitals/bands near the Fermi level. On the other hand the 
same mechanism can have a slightly detrimental effect at 
larger 3 J — U, where the ideal inter-orbital pairing condi- 
tions along the diagonals are weakened by the additional 
hybridization. One may also wonder if the Bogoliubov QP 
dispersions have anisotropic gaps. The resulting QP bands 
are shown in fig. [4] and are fully gapped with a fourfold 
symmetric gap modulation in k space, even though the 
gap minima are tiny. 

Note that the present superconducting state does not 
break TRS. The magnetic response is a combination of 
paramagnetic (spin-triplet) and spin-singlet behaviours, 
with a slightly larger out-of-plane than in-plane total 
magnetic susceptibility as shown in fig. [5j where M = 
(Li) + 2 (Si) is the total magnetization including orbital 
and spin contributions and Hb = B • 5^(1^ + 2S^) cou- 



ples the orbital and spin degrees of freedom to the exter- 
nal field B. Both orbital and spin expectation values are 
finite with roughly equal contribution to the total mag- 
netization. For comparison, the normal state magnetiza- 
tions are also shown in fig. [5] and are larger than in the 
superconducting state, as expected for a combination of 
spin-singlet and -triplet pairing in the presence of SO in- 
teraction. In particular, note that the spin magnetization 
changes drastically in the superconducting state with in- 
creasing A. In general, the magnitude of the d- vectors, and 
thus the magnetic response, can be modified by changing 
the size of the FS sheets. For instance a larger overlap 
between yz and xy dominated portions of the FS would 
enhance d yz / xy compared to d yz / xz and d xz / xy . The spin 
susceptibility then would be mostly dominated by d yz / xy , 
a situation which may be facilitated by applying uniaxial 
pressure. 

Discussion and summary. — Given that we based 
our MF study on the Sr2Ru04 compound to illustrate the 
effect of SO interaction on pairing, let us comment on the 
compatibility and the limitations of our results with what 
is known about the superconducting state in Sr2Ru04 [5j 
E]. Based on the QP gap variation along the FS sheets, 
one expects that this modulation may also be reflected 
in orientation sensitive specific heat measurements. Such 
magnetic field dependent specific heat measurements on 
Sr2Ru04 have indeed been carried out [31,32, but the 
interpretation of the experimental results is controversial, 
making a link to our QP dispersion difficult. However, due 
to the nature of inter-band pairing, the superconducting 
state presented here is sensitive to any kind of impurities 
associated with inter-band scattering, which is consistent 
with the phenomena observed in Sr2Ru04. 

Our result on the magnetization indicates that the spin- 
susceptibility is finite and different for in-plane and out- 
of-plane magnetic field orientations in both the normal 
and the superconducting state, as reported on Sr2Ru04. 
Yet below T c the in-plane and out-of plane susceptibilities 
decrease, which is in contrast to NMR Knight shift mea- 
surements [33j[3l], which revealed that a change in the 
spin-response across T c is absent for any field orientation. 
This behaviour differs also from the response expected of a 
chiralp + zp superconductor, where the spin-susceptibility 
decreases for field directions perpendicular to the a-b plane 
but remains constant for parallel orientations. While the 
amount of change in the present model depends sensitively 
on the SO interaction strength, as shown in fig. [5j the 
question also arises as to how orbital and spin contribu- 
tions were separated to obtain the Knight shift data when 
SO interaction is significant. Besides this, we note that 
the magnetic field effect on vortices will be highly non- 
trivial as well, as it involves competition between various 
types of vortices including half-quantum vortices [35| f36] 
in the presence of moderate SO coupling. 

Finally, the lack of TRS breaking is compatible with the 
absence of chiral supercurrents as observed in scanning 
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Hall probe and scanning SQUID measurements [371 138]. 
However, this contrasts with another proposal that the 
chiral states due to p + ip pairing on a and f3 bands cancel 
each other leading to a topologically trivial superconduc- 
tor [27]. It also contradicts Kerr rotation and /iSR mea- 
surements which have been interpreted in favour of TRS 
breaking [39,40 . The issue as to whether TRS is broken 
or not is not yet resolved in the experimental commu- 
nity. While the current study supports a non-TRS break- 
ing state, it can be modified by going beyond local interac- 
tions. A natural extension would be to include the effect of 
further neighbour ferromagnetic interactions such as those 
discussed by Ng and Sigrist [18], which could lead to a 
small admixture of odd parity pairing with broken TRS 
in addition to the pairing found here and which may be re- 
sponsible for the broken TRS signatures found in fiSK and 
Kerr experiments 39.140]. Another possibility is a finite- 
momentum pairing state such as a FFLO (Fulde-Ferrell- 
Larkin-Ovchinnikov) state JJTJ22]. It is plausible that a 
FFLO state between different bands can be stabilized over 
the inter-band pseudospin-triplet pairing. These studies, 
and more definite predictions for S^RuC^ or other specific 
materials, however, go beyond the scope of the current 9 
complex order parameter minimization and require more 
detailed work. 

In summary, we studied the combined effect of Hund's 
and SO coupling on t2 g orbital systems. Three orbital- 
singlet spin-triplet pairings were found to form an or- 
thogonal d-vector triad. A linear coupling between even- 
parity inter-orbital spin-triplet and even-parity intra- 
orbital spin-singlet pairings was allowed due to SO inter- 
action, determining the orientation of the three d- vectors 
and giving rise to a relative phase difference of tt/2 be- 
tween spin-singlet and spin-triplet order parameters. We 
also showed that inter-orbital spin-triplet pairing in the 
orbital basis corresponds to ever-parity inter- and intra- 
band pairing in the Bloch band basis, and discussed how 
the pairing strength varies within the Bloch bands. We 
further found that SO coupling assists Hund's coupling 
driven pairing, which generally leads to an anisotropic QP 
gap and an orbital dependent magnetic response. 
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